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ABSTRACT 


A fourth order Runge-Kutta method is applied to a system of linear 
differential equations. The solution over one time step is obtained 
explicitly as a function of initial conditions and the forcing function 
evaluated at the initial time, the final time, and the midpoint. 
Successive applications of this result can be employed to obtain an 
explicit solution for a number of time points which can later be used 
to rapidly determine solutions for a large number of different sets of 
initial conditions and forcing functions. 



NASA - GEORGE C 


MARSHALL SPACE FLIGHT CENTER 



NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 


* 


TECHNICAL MEMORANDUM X-53495 


EXPLICIT RDNGE-KUTTA INTEGRATION 

By 

Lyle R. Dickey 


TECHNICAL AND SCIENTIFIC STAFF 
AERO-ASTRODYNAMICS LABORATORY 
RESEARCH AND DEVELOPMENT OPERATIONS 



TECHNICAL MEMORANDUM X-53495 


EXPLICIT RUNGE-KUTTA INTEGRATION 
SUMMARY 


A fourth order Runge-Kutta method is applied to a system of linear 
differential equations. The solution over one time step is obtained 
explicitly as a function of initial conditions and the forcing function 
evaluated at the initial time, the final time, and the midpoint. 
Successive applications of this result can be employed to obtain an 
explicit solution for a number of time points which can later be used 
to rapidly determine solutions for a large number of different sets of 
initial conditions and forcing functions. 


EXPLICIT RUNGE-KUTTA 


Quite often it is necessary to solve the same system of differential 
equations for a large number of different initial conditions and forcing 
functions. In the event that the differential equations are linear, 
considerable computer time can be saved by taking advantage of the fact 
that the Runge-Kutta method of integration is a linear operation. With 
a little additional effort, the operation can be carried out explicitly 
resulting in a solution which is a linear combination of values of the 
forcing function evaluated at discrete points and the initial conditions. 
Having obtained the solution in this form, it is a simple matter for a 
computer to rapidly evaluate the solution for an extremely large number 
of different conditions. If this can be accomplished over one time step, 
it is only a matter of bookkeeping to repeat the operation over a number 
of time steps. With this in mind, the following system of linear differ- 
ential equations will be solved explicitly over one such time step by 
fourth order Runge-Kutta: 


X = F(X, t) = A(t) X + H(t) . 


( 1 ) 


The solution at a point t x can be obtained from the initial conditions 
at t c and the forcing function, H(t), evaluated at t Q , t x and the midpoint 
t by the following method: 



Xi = Xo + i (Ko + 2K i + 2K 2 + K 3 ), 


where 


Kq = hFO^, t Q ) 


K x = hF (xq + Y , t J 


K 2 = hF (x<) + y , 


K 3 = hF(Xo + K 2 , t x ) 

Xo - X(t Q ) 

X x = X(t x ) 
h = t x - t Q . 


( 2 ) 


Since the vector equation (1) is linear and the operations described 
in equations (2) are also linear, the solution is of the following form: 


X x = UXq + C q Ho + CH + Cx Hi , 


(3) 


where 

Ho - H(t c ) 

H = H(t) 

Hx = H(tx). 
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Successive substitution of equation (1) into the expressions 
described in equations (2) yield the following equations for U, C Q , C, 
and Cx* 


where 


U = I + t + 4A + A x ) + - 7 - (A Aq + A 2 + A x A) 


+ (A 2 Ao + A x A 2 ) + ~ <M A 2 Ao), 


C a, + tfi + £A^ilA, A 2 
c o 6 - 1 6 12 24 x 


2h 


h 3 


C = Y I + ^ (A + A x ) + J2 A l A 


C 1 - I I. 


Ao = A(t 0 ) 
A = A(t) 

Ax 35 A( tx ) • 


If the function H(t) is in tabular form and available at time incre- 
ments which differ from h, some sort of interpolation method is required. 
However, since most of the useful interpolation formulas are linear with 
respect to the dependent variable, H c , H and H x can be expressed as a 
linear combination of the values of H(t) at the points for which they are 
tabulated. Although the bookkeeping may become somewhat tedious when 
the solution is desired at a large number of time points, the use of 
matrix notation and digital computers, both of which are well designed 
for such operations, makes this an extremely small point for consideration. 
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The explicit method outlined above is extremely useful when a large 
number of different initial conditions and forcing functions are to be 
used. In particular, if the statistics of these inputs are available, 
the mean and standard deviation of the output can be evaluated in a 
straightforward manner and achieve results that would require an extremely 
large amount of computer time if the differential equations were to be 
solved over and over again for each individual case. 
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